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ABSTRACT 



We present observations and an interpretative model of the dust environment 
of Main-Belt Comet P/2010 F5 (Gibbs). The narrow dust trails observed can be 
interpreted unequivocally as an impulsive event that took place around 2011 July 
1st with an uncertainty of ±10 days, and a duration of less than a day, possibly 
of the order of a few hours. The best Monte Carlo dust model fits to the observed 
trail brightness imply ejection velocities in the range 8-10 cm s _1 for particle sizes 
between 30 cm and 130 /im. This weak dependence of velocity on size contrasts 
with that expected from ice sublimation, and agrees with that found recently for 
(596) Scheila, a likely impacted asteroid. The particles seen in the trail are found 
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to follow a power-law size distribution of index Ri-3.7. Assuming that the slowest 
particles were ejected at the escape velocity of the nucleus, its size is constrained 
to about 200-300 m in diameter. The total ejected dust mass is > 5 x 10 8 kg, 
with represents approximately 4 to 20% of the nucleus mass. 

Subject headings: Minor planets, asteroids: individual (P/2012 F5 (Gibbs) - 
Methods: numerical 



1. Introduction 



Main-Belt Comet P/2012 F5 (Gibbs) (hereafter P/Gibbs) was dis covered in the course 
of the Mt. Lemmon Survey on UT 2012 Mar. 22.89 JGibbs et alJboij ). The object showed 
a narrow dust tail 7' in length, and was identified as a Main-Belt Comet (MBC) owing to 
its orbital parameters. MBCs have dynamical properties of asteroids (i.e., Tisserand param- 
eters respect to Jupiter larger than 3), but physical characteristics of comets (gas and/or 
dust emission). The total members of this new class of Small Solar System Objects amount 
to nine (the 9th becoming P/Gibbs), so that the statistics is still poor so as to identify 
th eir globa l phys ical p roperties or d ynamical history. For reviews on those objects, we refer 



to iBertinil (j201l[ ) and iJewittl ( 120121 ) . Most of the MBCs have been shown to be dyna mi- 
cally stable on timesc ales of 100 Myr or longer ( Ipatov fc Hahn 1999 : Hsieh et al. 2012a b cl; 
Stevenson et all2012 ) suggesting that they are native me mbers of the Main Asteroid Belt and 
not captured objects from elsewhere (IHsieh et al.ll2009l ). This is strongly supported by dy- 
namical and spectroscopic arguments that show that some of them belong to well established 
aste roid collisional famili es and that their spectra is different to that of "normal" comets 



(e.g. iLicandro et al.ll201ll ). Howeve r, there are some members as 238P/Re ad and P/2008 Rl 
that are stable for 20-30 Myr only ( Haghighipour 2009 : Jewitt et al. 12009). Regarding thei r 
activity, some are clearly recurrent, as 133P/Elst-Pizarro and 238P (IHsieh et al.ll2010l . 1201 ll ). 
while some others like (596) Scheila appear to have ejected dust via an imp u lsive event that 



20 111 : iMoreno et al.l 12011 at llshiguro et al. 



migh t be associated to a collision (I Jewitt et al.l 12011c lYang &: Hsiehl 1201 it iBodewits et al. 



20111 ). On the other hand, P/2010 R2 (La Sagra) 



and 2006 VW139 appe ared to be continuous l y active during periods of at least 200 d ays and 
100 days, respectively ( Moreno et al. 2011b : Licandro et al. 2012 : Hsieh et al. 2012bl ). 



Observations and analysis of P/Gibbs have been recently shown by IStevenson et al. 



( 120121 ). who report on the impulsive character of the emission of the dust, the outburst 
date, the mass of the ejecta, and establish limits to the nucleus size. In this paper we 
report images in the red spectral domain acquired with instrumentation attached to the 
10.4-m Gran Telescopio Canarias (GTC) of P/Gibbs at two different dates, and perform an 
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interpretation of the observed trail brightness using a forward Monte Carlo tail model, in 
order to provide estimates of the dust emission times, the ejected mass, the particle ejection 
velocities, and their size distribution function. We also set limits to the nucleus s ize based 
on the escape velocity, and compare our results to those by lStevenson et al.l ( 120121 ). 



2. Observations and data reduction 

CCD images of P/Gibbs were collected under photometric conditions on the nights of 
18 May and 8 June 2012, using a Sloan r' filter in the Optical Syst em for Image and Low 



Reso lution Integrated Spectroscopy (OSIRIS) camera-spectrograph ( ICepa et al.ll2000t ICepa 



20101 ) at the GTC. The OSIRIS instrument consists of two Marconi CCD detectors, each 
with 2048x4096 pixels and a total unvignetted field of view of 7.8' x 7.8'. The plate scale was 
0.127 "/px, but we used a binning of 2x2 pixels in order to improve the signal to noise ratio, 
so that the spatial resolution of the images is 485 km px -1 and 540 km px" 1 on May 18.9 
and June 8.9, respectively. The images were bias subtracted and flat-field corrected using 
standard procedures. A total of 18 images were acquired each night and calibrated using 
standard stars. The images were converted to solar disk intensity units appropriate for the 
analysis in terms of dust tail models, and a median stack image was obtained for analysis for 
each date (see Figure 1). As a result of both the flux calibration and the median stacking 
procedure, we estimated a total flux uncertainty in the com bined images of 0.3 mag. The 



images were finally rotated to the (N,M) coordinate system (IFinson fc Probsteinlll968l ). 



The Model 



We have performed an analysis of the two images shown in Figure 1 by a direct Monte 
Carlo dust tail model. This model was used to characterize the dust environments of various 



comets and MBCs (e.g.. lMorenoll2009uMoreno et al.ll2010l . l2011bl ). being applied in particular 
to comet 67P/Churyumov-Gerasimenko, th e target of Rosett a Mission to arrive to the comet 
in 2014 (the so-called Granada model, see iFulle et al.l 120101 ) . The code is used to compute 
the trajectory of a large number of grains ejected from a cometary or asteroidal surface. We 
assume that the only governing forces on those dust particles are the solar gravity and the 
radiation pressure. Thus, the gravity of the object itself is neglected, which constitutes a 
good approximation for small-sized nuclei. In the case of cometary activity, we consider that 
the particles are accelerated by gas drag from ice sublimation to their terminal velocities, 
which are the input ejection velocities considered in the model, that is applicable for any 
other ejection mechanism. Once ejected, the particles describe a Keplerian trajectory around 
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the Sun, whose orbital elements are computed from the terminal velocity and the ratio of 
the force exert ed by the solar radiation pressure and the solar gravity (the /3 parameter, 



sec 



Fulldll989l ). This parameter can be expressed as (3 = C pr Q pr /(2pr), where C pr =1.19x 



10 kg m , Q pr is the radiation presure coefficient, and p is the particle density. For each 
observation date, the trajectories of a large number of dust particles are computed, and then 
their positions on the (N, M) plane are calculated. Finally, their contribution to the tail 
brightness is computed. 

Before applying the model, we performed a preliminary analysis of the images in terms 



of a synchrone map (IFinson fc Probsteinlll968l ). which provides us with an approximate idea 



on the time interval at which the particles were ejected from the nucleus related to the 
observation date. Figure 2 shows the two images as well as the corresponding synchrones at 
times of ±100, ±20, ±10 and days relative to the 2011 July 3rd synchrone, the one which is 
approximately best aligned with the narrow trail in each image. The smallest sized particles 
in the 2011 July 3rd synchrones have r ~130 /iin and r ~180 /xm for the two observing 
dates, respectively. These diagrams indicate that all the dust particles having r >130 /zm 
were ejected during a very short time interval around that date, because otherwise the 
trails would have been wider towards the lowermost portions of the trails, even if they were 
ejected with zero velocity respect to the nucleus. Therefore, an impulsive event must be the 
responsible of the ejection of dust, so that an outburst, a collision with another body, or a 
rotational disruption could be in principle invoked among the responsible mechanisms. The 
precise event date, and its duration, is to be determined on the basis of the Monte Carlo 
dust tail analysis (see section 4). 

An important aspect of the observations and the modeling is the location of the asteroid 
nucleus. We will work under the assumption that the object nucleus is immersed in the dust 
cloud, and that the optocenter of this cloud corresponds to the actual location of the nucleus. 
We further assume that the observed brightness is dominated by the dust, the contribution 
of the nucleus being negligible. This hypothesis would be tested on the basis of the nucleus 
size derived from arguments based on the escape velocity. 

A number of simplifying assumptions on the physical parameters must be made in 
order to make the problem tractable. Thus, we assume that the trail is composed of spher- 
ical particles of carbonaceo us composition, having a refractive index at red wavelengths of 



m=1.88+0.71i (jEdohl Il983l ) . which implies a geometric albedo of p„=0.036, and a pr essure 



radiation coefficient of Q pr ~ 1 for particles of radius r >1 /im (IMoreno et al.l 120121 . their 
Figure 5). We further assume that the particles have a density of 1000 kg m -3 . The particle 
size distribution is assumed to follow a power- law of index a, and the terminal velocities are 
described as a function of the (3 parameter as v(/3) = vq/3 1 . This relationship is generally 
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accepted for the termin al velocities of comet dust, and also f or fragments ejected from col- 
lision experiments (e.g., lGiblinlll998l ; lOnose &: Fujiwaral 12004 ) . The total ejected dust mass 
(M e ) and the maximum size of the ejected particles (r max ) must be also specified, as well as 
the starting time (t s ) and the duration of the event (At). In summary, there are a total of 
seven adjustable parameters to fit the observations (a, vq, 7, M e , r max , t s , and At). 



4. Results 

We attempted to fit the observed dust trail brightness by minimizing the function a = 
&\ + V2, where the subscripts 1 and 2 correspond to the images obtained in 2012 May 
18.9 and 2012 June 8.9, respectively, and <7j = a/ (^JJobJj) ~ I fit{i)) 2 / N (i)) , where I b s (i) 
are the observed trail intensities and Ifu{i) are the fitted intensities, the summation being 
extended in principle to all the image pixels N(i). Since there are some regions in the images 
that are strongly contaminated by bright field stars, we restricted the summation to pixels 
outside those regions, and located mostly along the trails. The minimization procedu re was 



performed by the multidimensional downhill s implex algor i thm (INelder &: Meadlll965l ). using 



the FORTRAN implementation described in iPress et al.l (119921 ) . Each of those parameters 



influence the derived trail brightness in different ways. Thus, Vq controls the width of the 
trail, 7 influences the variation of the width of the trail along it, a constraints the slope of 
the brightness along the trail, the total ejected mass influences the overall brightness, and 
the maximum particle size controls the brightness mainly at the head of the trail, so that 
if f^max is set to a small value (say, smaller than several centimeters) the peak of brightness 
is displaced significantly trailward and the fit becomes impossible. Regarding the minimum 
particle size, it is not a free parameter, as it is found by the intersection of the synchrone 
that is best aligned with the trails and the N axis, which corresponds in our images to ~130 
fim and ~180 fim for 2012 May 18.9 and 2012 June 8.9, respectively. 

Since the downhill simplex method searches for a local minimum of the function in the 
parameter space, we performed several runs by varying the starting simplex. We found that 
all runs tend to converge to close local minima, giving similar values of a. The deepest 
global minimum was then taken as the best fit, which has <r=4.9xl0 -15 solar disk units (see 
Table 1, and Figure 3). Unfortunately, this technique does not provide any estimate of the 
uncertainty in the derived parameters. To estimate those uncertainties, we determined the 
errors by perturbing the observed image intensities by the flux uncertainties and then finding 
a new best fit. In addition, we verified that outside the error limits for each one of the fitted 
parameters displayed in Table 1, no satisfying solutions were found. To do that, we tried to 
fit the model with values outside the parameter bounds specified in the table, and could not 
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find satisfying solutions in any case. 



An important result is that the derived ejection velocities are almost independent on 
size (7=0.041002)5 ranging from about 10 to 8 cm s _1 for particles between the lower and 
upper size limits (130 fim and 28 cm). This flat dependence of ejection velocities on size 
have been previously inferre d in our analysis of t he outburst of the asteroid (596) Scheila, 
for which we found 7=0.05 ( jMoreno et al.ll2011al ) . This dependence is markedly different 
to that expected from gas drag by ice sublimation processes as occur in most comets, with 
typical values of 7 ~0-5, and could be associated to a collision event as was suggested for 



Scheila ( jMoreno et al.l l2011aJ ) . although it cannot be confirmed. On the other hand, if we 
equate the velocity of the slowest moving particles (v=8 cm s _1 for r=28 cm particles) to 
the escape velocity, we obtain a nucleus radius of R n =107 m to i? n =152 m, considering 
bulk densities of 1000 to 500 kg m -3 . At this point, it is interesting to note that computer 
simulations at the catastrophic disruption threshold reveal that ejecta velocities depend on 
target size and that ejection velocities of the order of those obtained here (~ 10 cm s~ l ) 



could be compatible with both porous and non-porous targets of ~100 m radius ( jJutzi et al. 



2010). Our d erived target rad i us is a far more stringent constraint than that of R n <2.1 
km derived by IStevenson et al.l ( 120121 ) based on the non-detection of the nucleus on images 
taken by the Wide-Field I n frared Survey Explorer (WISE). On the other hand, using the 
formalism by iBowell et al. rll989h . and assuming a slope parameter of G=0.15, and a bulk 
geometric albedo of p„=0.15, the apparent magnitude of a i? n =125 m nucleus at the same 
phase angle and geocentric and heliocentric distances than P/Gibbs would be just above 
m„=26, which represents a negligible contribution to the brightness compared with that at 
the optocenter of the trails. This is compatible with our initial hypothesis on that the dust 
cloud brightness dominates the trail optocenter. An object of such a small size might have a 
very small rotation period that could be under the critical rotation pe riod for fract ure so that 
a rotational disruption could also be argued as the cause of the event ( ]Jewittll2012l ). provided 
it produces a sudden release of material. A small-sized nucleus was also derived for P/2010 
A2 (LINEAR), the innermost MBC discovered so far, whose tail has been also reported by 



many authors as the resu lt of an impulsive event (e.g. IJewitt et al.l |2010| ; ISnodgrass et al. 



2010 



f lMoreno et al.ll2010h 



Hainaut et al.l 120121 ). and not as the result of a sustained activity as we suggested 



The size distribution power index, which essentially controls the slope of the trail bright- 
ness, is -3.7±0.1. The maximum ejected particle size derived is r=28±10 cm (Table 1). A 
maximum size smaller than the lower bound would imply a significant displacement of the 
maximum of brightness trailward. Conversely, a size larger than the upper bound would 
make the slope of the synthetic trail much steeper than observed. 
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The minimum e jected dust mass is (5±2)xl0 8 kg, which is near the value reported by 
Hainaut et al.l (20121) for P/2010 A 2 (8xl0 8 kg). This is an order of magnitude higher than 
reported by IStevenson et al.l (120121 ). although their result is just based on a measure of the 
total geometric cross section of the trail and assuming a mean size for the ejected particles, 
not a size distribution. This total mass is in fact a lower limit, since we cannot precise 
neither the amount of small (r <100 /zm) particles released, nor the component of particles 
traveling at high speed that could result as consequence of e.g. a collision, if this were the 
case. 



The best-fit event date (2011 July 1st) agrees with that derived by IStevenson et al. 



( 120121 ) within the error limits (they derived 2011 July 7th as the event date). The duration 
of the event is constrained to be of less than a day. In the At <1 day space, it is always 
possible to find a set of fitting parameters giving similar values of a. The fits shown in Figure 
3 correspond to an event duration of 0.13 day, i.e., just above 3 hours. An impulsive event 
is, then, clearly implied. 



5. Conclusions 

From the Monte Carlo dust tail modeling of the observations of Main-Belt Comet P/2012 
F5 (Gibbs) we can derive the following conclusions: 

1) The direct study of the MBC images in term of synchrone analysis imply an impulsive 
event as the cause of the observed trails. We predict the event to have occurred on 2011 
July 1st, wi th an accuracy o f ±10 days. In this respect, we agree with the results recently 



reported by IStevenson et al.l (120121 ) for this MBC, not only in the nature, but also in the 
event date within the errors (they report on a impulsive event on 2011 July 7th with an 
uncertainty of ±20 days). The duration of the event is constrained to be less than 1 day, 
possibly less than a few hours. 

2) The nature of the impulsive event is impossible to determine with the sole information 
of the physical parameters derived from this analysis. Among the likely causes, an outburst, 
a collision with another body, or a rotational disruption could be invoked. Activity related 
to ice sublimation seems unlikely on the basis of the dependence of the particle ejection 
velocitie s on size, which inter estingly turns out to be very similar to that we found for (596) 



Scheila ( iMoreno et al.l 1201 laf ). a very likely impacted asteroid. 



3) The total dust mass released is >5xl0 8 kg. The ejected particles are distributed in 
size following a power-law of index -3.7±0.1. The maximum particle size ejected is about 
30 cm, with a velocity of ~8 cm s _1 . Adopting this value as the escape velocity, the size of 
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the nucleus is constrained to about 100-150 m in radius, for bulk body densities in the range 
1000-500 kg m" 3 . 

This article is based on observations made with the Gran Telescopio Canarias (GTC), 
installed in the Spanish Observatorio del Roque de los Muchachos of the Instituto de As- 
trofisica de Canarias, in the island of La Palma. 
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Fig. 1. — Median stack images of P/2012 F5 (Gibbs) obtained with the OSIRIS instrument 
of the 10.4m Gran Telescopio Canarias through a Sloan r' filter, on UT 2012 May 18.9 (upper 
panel) and June 8.9 (lower panel). The directions of the velocity vector, the Sun, and the 
astronomical North and East are indicated. 
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Fig. 2.— Panel(a) displays the P/2012 F5 (Gibbs) image obtained in 2012 May 18.9. Panel 
(b) shows the same image but with seven synchrones overlaid, corresponding, in clockwise 
order, to -100, -20, -10 (dotted lines), (solid line), +10, +20, and +100 (dashed lines) 
days respect to the 2011 July 3 synchrone. Panels (c) and (d) shows the same as (a) and 
(b), respectively, but for the 2012 June 8.9 image. The scale is 485 km px" 1 , and 540 km 
px _1 , in the upper and lower panels, respectively. Note that the x and y axes are scaled 
independently to facilitate comparison between synchrones. 
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Fig. 3. — Panels (a) and (b) correspond to the observation and best fitted image for the 
observation on 2012 May 18th. Panel (c) shows observed and modeled scans along the trail 
of those images. Panels (d), (e), and (f) give the same as (a), (b), and (c), respectively, 
but for the 2012 June 8th image. The spatial scales are 485 km px -1 , and 540 km px -1 , 
respectively. Note that the images are rotated by 90 degree counterclockwise so that their x 
and y axes correspond to the photographic M and N axes. The segments labeled as "BS" 
pertain to regions of field star contamination in the trails. 
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Table 1. The best-fit parameters of the model. 
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